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Comparison of the full-potential and frozen-core approximation approaches 
to density-functional calculations of surfaces 
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We scrutinize the accuracy of the pseudopotential approximation in density-functional theory 
(DFT) calculations of surfaces by systematically comparing to results obtained within a full-potential 
setup. As model system we choose the CO oxidation at a RuO2(110) surface and focus in particular 
on the adsorbate binding energies and reaction barriers as target quantities for the comparison. 
Rather surprisingly, the major reason for discrepancy does not result from the neglected semi-core 
state relaxation in the frozen-core approximation, but from an inadequate description of the local 
part of the Ru pseudopotential, responsible for the scattering of / like waves. Tiny, seemingly 
irrelevant, imprecisions appearing in these properties can have a noticeable influence on the surface 
energetics. At least for the present example, we obtain excellent agreement between both approaches, 
if the pseudopotential describes these scattering properties accurately. 

PACS numbers: 68.43.Bc, 71.15.Dx, 82.65.+r 
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I. INTRODUCTION 

Density-functional theory (DFT) calculations have be- 
come an important tool in present-day materials science. 
This does not mean, however, that such calculations are 
routine or should be handled as a black-box. On the 
contrary, the large scatter in published DFT numbers 
demonstrates the crucial importance of how the calcula- 
tions are technically carried out, and to scrutinize which 
approximations are employed to make computationally 
most intensive applications like large-scale surface stud- 
ies feasible. The pseudopotential approach based on the 
discrimination between core and valence electrons is one 
such approximation to significantly reduce the compu- 
tational burden 1 . It considers the chemically inert core 
electrons together with the nuclei as rigid non-polarizablc 
ionic cores, so that only the valence electrons have to 
be dealt with explicitly. Continuing efforts to improve 
the accuracy of this approximation, while maintaining 
or even further reducing the computational cost, have 
first led to the development of ultra-soft pseudopoten- 
tials (USPPs) 2 and ultimately to the formulation of the 
projector augmented wave (PAW) method^. 

Over the years, the pseudopotential approximation has 
proven its accuracy and utility for calculations of bulk 
properties for a variety of materials with widely different 
bonding character. In fact, the approximation has been 
so successful that nowadays database pseudopotentials, 
provided through the large program packages, are rou- 
tinely used. If the quality and suitability of these pseu- 
dopotentials for a specific study is then checked at all, 
it is normally deemed sufficient to perform a few com- 
putationally inexpensive tests. Typically this could in- 
clude looking at atomic scattering properties or to com- 



pare the results for standard bulk properties like lattice 
constants against those obtained in full-potential calcula- 
tions or those measured experimentally. Obviously, such 
tests become less reliable, the more the target quanti- 
ties of interest deviate from those checked. Furthermore, 
the comparison to other theoretical studies becomes ob- 
scured by the other approximations and uncertainties in 
practical DFT computations like differently converged 
basis sets, while comparison to experiment does not ac- 
count for the error introduced by the approximate DFT 
exchange-correlation (XC) functional. 

This problem applies prominently to calculations at 
surfaces or for surface reactions. Due to the large com- 
putational costs of such calculations, the pseudopotential 
approximation has rarely been tested with sufficient care. 
Published DFT numbers often show such an unsatisfy- 
ing scatter that even the use of different XC functionals 
is sometimes less decisive than the use of other approx- 
imations as, for example, pseudopotentials, finite basis 
sets or k-point meshes. In the present study we there- 
fore set out to carefully check on the accuracy of the 
pseudopotential approximation for surface calculations. 
This is done by systematically comparing to results ob- 
tained within the full-potential linear augmented plane 
wave (LAPW)i and augmented plane wave plus local 
orbitals (APW+lo) method^. As target quantities for 
the comparison, we focus on adsorbate binding energies 
and reaction barriers as typically of interest in hetero- 
geneous catalytic systems. Since transition metals are 
frequently used in such applications, we also analyze the 
role of so-called semi-core states, i.e. states that lie en- 
ergetically intermediate between "true" core and valence 
states. Since the wavefunctions of semi-core states expe- 
rience some (though small) overlap with their counter- 
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parts of neighboring atoms, they need special attention 
in both full-potential and frozen-core approaches. In the 
LAPW/APW+lo method they are treated by a local or- 
bital basis se'ba, while in the PAW method they can either 
be kept frozen in the core or included in the valence states 
at a slightly higher computational cosA 

As a model example for our test we consider CO oxi- 
dation at the RuO2(H0) surface. This system has re- 
cently attracted a lot of attention as a highly active 
model catalyst&L&ii, in particular after extensive experi- 
mental and theoretical work had shown that an epitaxial 
RuO2(110) film forms at the also much studied Ru(0001) 
model catalyst surface at realistic O2 pressure a 10 ! 11 ! 12 ! 13 . 
Rather than by this interesting physics, our choice is, 
however, motivated by the large number of performed 
DFT studieo 7 ! 14 ! 15 ! 16 ! 17 ! 18 ! 19 ! 20 . Symptomatic for many 
surface systems the published adsorbate binding energies 
show a significant scatter of up to 1 eV, even though they 
were computed with the same XC functionali^i^iSiiLi^. 
Our analysis, reported below, reveals that a major source 
for discrepancies is caused by an inadequate description 
of the Ru 4/ scattering properties, and, if special care is 
taken in the pseudopotential creation, the neglect of core 
and semi-core relaxation is found to play only a minor 
role. 

We believe that this often overlooked sensitivity of the 
pseudopotential approximation is of wider importance 
than just for the RuO2(110) system studied here. Apart 
from insufficiently converged basis sets, such hitherto not 
much addressed deficiencies in the employed pseudopo- 
tentials can be a crucial factor behind the often substan- 
tially differing results reported in the DFT literature of 
this and other surface systems. 



II. COMPUTATIONAL DETAILS 

Our study focuses on comparing full-potential results 
to those obtained within the pseudopotential approxi- 
mation. For both methods we choose widely used ap- 
proaches, namely the LAPW/APW+lo scheme^ as im- 
plemented in the WIEN2k cod o 21 i 22 for the prior and the 
PAW method 3 as implemented in the Vienna ab-initio 
simulation package (VASP) 2 ^ 2 ^ for the latter. To make 
this comparison meaningful requires to make the two 
kinds of calculations as similar as possible. We there- 
fore use exactly the same geometric structures, exactly 
the same k-points and exactly the same XC functional. 
If care is then taken to have convergence with respect 
to the remaining basis set parameters, all remaining dif- 
ferences must be entirely due to the different treatment 
of the core and semi-core electrons in the two methods. 
Specifically, we increased the employed basis sets until 
the targeted surface quantities, i.e. adsorbate binding 
energies and reaction barriers, were absolutely converged 
to within ±10meV. 

Here, we first detail those computational parameters 
that are essentially identical in both approaches, and 



then describe in the next two sub-sections those param- 
eters that are treated differently. The RuO2(H0) sur- 
face was modeled in the supercell approach by symmet- 
ric slabs, consisting of three 0-(Ru202)-0 trilayers and 
separated by approximately 11 A vacuum. The O atoms 
and CO molecules are adsorbed symmetrically on both 
sides of the slab. All atomic positions in the outermost 
trilayers were fully relaxed in the course of previous full- 
potential calculations^, and they were kept fixed as in- 
put positions for all calculations presented in this work. 
Both the LAPW/APW+lo and the PAW calculations 
were thus performed using exactly the same structural 
parameters. Adsorption of O and CO is considered in 
(1 x 1) surface unit-cells, while all reaction barriers are 
obtained using larger (1 x 2) cells as shown in Fig. Elbe- 
low. Dense (5 x 10 x 1) and (4x4x1) Monkhorst-Pack 
special k-point meshes were used to sample the Brillouin 
zones of these cells, respectively. For the fractional oc- 
cupation numbers a temperature smearing of 0.2 eV was 
applied in LAPW/APW+lo (with the corresponding ex- 
trapolation to T -> OK 2 ^), while in the PAW calcu- 
lations the Methfessel-Paxton method^! with a width of 
0.2 eV was used. All calculations utilized the exchange- 
correlation functional based on the PBE version of the 
generalized gradient approximation 2 ^. 



A. LAPW/APW+lo method 

In order to achieve an accurate description without 
any shape approximation for the potential the LAPW 
and APW+lo methods use specially adapted basis sets. 
Inside non-overlapping, so-called muffin tin (MT) spheres 
centered around the atomic sites, in the LAPW method 
linear combinations of solutions of the radial Schrodinger 
equation and their energy derivatives are employed to 
augment a plane wave basis set in the remaining, so- 
called interstitial region between the MT spheres^. In 
the more recent APW+lo methocU, this augmentation is 
done by only using linear combinations of solutions to the 
radial Schrodinger equation and providing extra flexibil- 
ity by adding local orbitals (lo) to the basis set. The lat- 
ter are radial functions confined to the MT sphere, which 
are also particularly useful to accurately describe semi- 
core states. As shown by Madsen and coworkers, using 
a mixed LAPW/APW+lo basis set for different angu- 
lar momentum values of radial functions centered at the 
same atom yields a particularly efficient description 29 , 
and this is what has been implemented into the WIEN2k 
cod o 21 i 22 . In the present work we correspondingly used 
an APW+lo basis for I = 0, 1, 2, while LAPW basis func- 
tions were used for all higher angular momenta up to 
^max — 12 for the wave functions and up to = 6 
for the potential. This means in particular, that local 
orbitals are set to specifically treat the Ru 4s and Ap 
semi-core states, as well as the O 2s orbital. 

The radii Rmt of the MT spheres around each ele- 
ment were chosen as follows: -Rmt = l-8bohr (0.95A), 
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R^ T = l.lbohr (0.58 A) and Rg T = l.Obohr (0.53 A). 
For this choice all remaining basis set parameters were 
optimized to achieve the aspired ±10meV absolute con- 
vergence of all later discussed binding energies and re- 
action barriers. Apart from the energy cutoff for the 
plane wave representation of the interstitial potential 
(i?pot x = 196 Ry), this concerns primarily the plane 
wave cutoff £™ f ax for the interstitial wave functions as 
the parameter most crucially determining the accuracy 
and computational demand. We were forced to increase 
this value up to 408 eV (30 Ry) for the chosen set of 
muffin tin radii. We note that such a cutoff would at 
present imply a prohibitive computational cost for most 
exploratory large-scale surface studies. This also explains 
the quantitative differences in the numbers reported be- 
low compared to preceding full-potential calculations for 
this systemi2*i2*2fi. In these studies, the DFT energetics 
served primarily as input to thermodynamic or statisti- 
cal mechanical approaches. The latter depended either 
on faster converging energetic differences^^ or were 
not much affected by modest errors in the underlying 
energetics^. The targeted convergence level was there- 
fore considerably lower, employing a 272 eV (20 Ry) cut- 
off in a pure LAPW basis. 



B. PAW method 

For our PAW study the frozen-core potentials were 
generated according to the procedure outlined in Ref. 0. 
The starting point for the investigations were the poten- 
tials provided originally within the VASP database. In 
these potentials, to which wc will from now on refer to as 
"PAW-std" , two partial waves were constructed for each 
angular momentum, and a pseudopotential truncated all- 
electron potential was used as local component. For oxy- 
gen and carbon the core radii for the s and p partial 
waves were r c (0) = 0.63 A and r c (C) = 0.80 A, respec- 
tively. For ruthenium the core radius for the 5s, bp and 
id partial waves was set to r c (Ru) = 1.30 A. The local 
pseudopotential was created by replacing the all-electron 
(AE) potential inside a sphere with a radius of 0.94 A 
by a single spherical Bessel-function, with the spherical 
Bessel-function chosen such that the resulting potential 
is continuously differentiable. 

To check on the effect of the various freedoms offered 
in the construction of PAW potentials, we also gener- 
ated several other pseudopotentials to compute all tar- 
get quantities used for the comparison. Most notable for 
the discussion below, these are two additional Ru PAW 
potentials, in which in one the 5/ states were chosen 
for the local potential ("PAW-/") and in the other the 
4p semi-core states were additionally treated as valence 
( "PAW-4p" ) . For the latter purpose one partial wave was 
constructed for the 4p and bp states, and simultaneously 
the radial cutoff for the 4p partial waves was reduced to 
1.10 A. Furthermore, we also employed very hard O and 
C potentials (0_h and C_h in the VASP PAW database) 



with a small core radius (r c = 0.58 A), which allows an 
accurate treatment of molecules like O2 and CO 30 , and 
corresponding reaction intermediates with similarly short 
bondlcngths. We add a "+" to the potential notation 
(e.g. "PAW-4p+"), when these hard pseudopotentials 
were used in the calculations. For the reference energy 
of the O2 molecule we always used the values calculated 
using these hard pseudopotentials. 

The energy cutoff required to achieve the aspired 
±10meV convergence of all energetic quantities used 
in the comparison is dictated by the oxygen potential. 
While routinely such calculations would be done with a 
400 eV (29.4 Ry) cutoff, we had to increase it to 600 eV 
(44.1 Ry) for the standard O PAW to match the tight 
convergence criteria imposed in this study. For the cal- 
culations involving the hard O and C PAWs a further 
increase up to 800 eV (58.8 Ry) was necessary. 



III. RESULTS AND DISCUSSION 

A. Standard pseudopotential tests 

We start our assessment of the frozen-core approxima- 
tion by subjecting the various PAW potentials to some of 
the computationally inexpensive tests often employed in 
the literature. Such necessary, but of course not sufficient 
tests address the logarithmic derivatives, ghost states, or 
transferability 3 -!^. Here, we focus on the logarithmic 
derivatives as a measure of the pseudopotential's ability 
to reproduce atomic scattering properties. For this one 
checks how well the logarithmic derivatives of the radial 
wave functions agree for the pseudo and the all-electron 
atom. This is evaluated as a function of the energy for 
the relevant s, p and d angular momenta at some diagnos- 
tic radius from the nucleus, say half a typical interatomic 
distance. FigureHshows this at a distance of 1.4 A away 
from the nucleus for the three Ru PAWs that we gener- 
ated and employed in our study. In the energy range of 
interest, i.e. the valence band regime, all three Ru PAWs 
pass this test, and describe the all-electron s, p and d 
scattering properties essentially exactly. There are only 
small deviations in the normally not much regarded / 
scattering properties of the standard PAW-std potential, 
which can be remedied by choosing an / pseudopotential 
as local component as in the PAW-/ potential, cf. the 
middle panel in Fig. ^ A remaining small discrepancy 
for the p scattering properties visible at energies clearly 
outside the valence band regime, i.e. below — 1.5Ry, can 
only be overcome by treating the Ap semi-core states as 
valence, as shown in the bottom panel in Fig. 2] (PAW- 
4p). 

Since tests like looking at the scattering properties 
assume spherical symmetry and neutral atoms, another 
class of frequently employed tests addresses non-spherical 
environments. For this, it had been noted that tests 
of the electronic hardness are particularly sensitiv e) 33 ! 34 . 
Another possibility is to cross check how well the pseu- 
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FIG. 1: Atomic scattering properties of the three Ru PAW 
pseudopotentials employed in the comparison (see text for la- 
beling and details of each potential): PAW-std (top panel), 
PAW-/ (middle panel) and PAW-4p (bottom panel). Shown 
are the logarithmic derivatives of the radial wavefunctions for 
different angular momenta for a spherical Ru atom, evaluated 
at a distance of r = 1.4 A from the nucleus. Solid lines corre- 
spond to the all-electron full-potential, and dotted lines to the 
PAW potential. The energy zero corresponds to the vacuum 
level. 



dopotentials reproduce some properties of simple test 
systems. Typically this involves bond lengths and bind- 
ing energies of molecules^, or lattice constants and for- 
mation energies of bulk crystals. Comparing against ex- 
perimental data is in this respect not of much use, since 
the approximations in the employed XC functional lead 
themselves to systematic deviations between theoreti- 
cal and experimental numbers, which need to be distin- 
guished from errors rooted in the pseudopotential itself. 



TABLE I: Comparison of structural and energetic properties 
of molecular and bulk test systems, obtained using the var- 
ious PAW potentials and within the LAPW/APW+lo full- 
potential method. Compiled are the bond length d and bind- 
ing energy E b of O2 and CO, as well as the lattice parameters 
a, c and the formation energy per formula unit Ef of bulk ru- 
tile Ru02- The binding energies are computed with respect 
to non-spherical spin-polarized atoms, and the formation en- 
ergies are calculated with respect to O2 and relaxed bulk hep 
Ru. 



o 2 


d(A) 




E b (eV) 


LAPW/APW+lo 


1.218 




6.21 


PAW-std 


1.232 




6.05 


PAW-std+ 


1.218 




6.22 


co 


d 




E b (eV) 


LAPW/APW+lo 


1.138 




11.65 


PAW-std 


1.143 




11.51 


PAW-std+ 


1.136 




11.65 


Bulk Ru0 2 


a (A) 


c(A) 


E f (eV) 


LAPW/APW+lo 


4.500 


3.117 


2.99 


PAW-std 


4.550 


3.129 


2.85 


PAW-/ 


4.525 


3.115 


3.00 


PAW-4p 


4.525 


3.115 


3.04 



In fact, it is frequently the "virtue" of a bad pseudopo- 
tential that it leads to excellent agreement with some 
experimental test data, despite the use of the approxi- 
mate XC functional. Instead, the comparison should be 
done with full-potential calculations using the same XC 
functional. Then, however, one needs to ensure that the 
comparison is not hampered by other approximations like 
insufficiently converged basis sets. 

Table ^ compiles such a comparison to the converged 
full-potential data from our LAPW/APW+lo computa- 
tions, both for the O2 and CO molecules and for the 
Ru02 bulk crystal. Again, the conclusion from such 
data would normally be that the O, C and all three Ru 
PAW potentials do a great job. The structural param- 
eters are all reproduced within 1 % error, and also the 
energetic parameters agree within a few percent. Al- 
though the Ru PAW-/ and PAW-4p potentials clearly 
outperform the standard PAW-std potential in the cal- 
culated Ru02 formation energy, the error per O-Ru bond 
amounts in the latter still only to 0.14/6 ~ 0.02 eV, con- 
sidering that every Ru atom in bulk rutile Ru02 is six- 
fold coordinated. Although experienced pseudopoten- 
tial practitioners know that errors in bulk tests multiply 
in low-symmetry environments, this deviation still looks 
rather small. The conclusion from these data and from 
the scattering property transferability test would there- 
fore most probably have been that all generated PAW 
potentials should be quite reliable, and could be used to 
obtain accurate surface energetics within the pseudopo- 
tential approximation at a much reduced computational 
cost compared to the full-potential approach. 
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cus site bridge site 




FIG. 2: Top view of the RuO 2 (110) surface showing the (1 x 1) 
and (1 x 2) unit-cells employed in the calculations, as well as 
the location of the two prominent adsorption sites, namely the 
bridge and the coordinatively unsaturated (cus) site. Shown 
is the surface termination exposing a (RU2O2) layer. Ru = 
light, large spheres, O = dark, medium spheres. Atoms lying 
in deeper layers have been whitened for clarity. 



B. O and CO binding at Ru0 2 (110) 

We now proceed to verify this positive assessment of 
the generated PAW potentials obtained from the preced- 
ing transferability tests by directly comparing surface en- 
ergetic quantities computed within the pseudopotential 
and the full-potential approaches. In a first step, this con- 
cerns the adsorption properties of simple probe molecules 
at the surface, which in the specific model system used 
translates into looking at the binding energies of O and 
CO at prominent adsorption sites at the RuO2(H0) sur- 
face. In the (110) direction, the bulk stacking sequence of 
rutile RuC-2 can be expressed as a series of 0-(Ru202)-0 
trilayersiiiS. Depending after which layer the sequence 
is truncated, three different (lxl) surface terminations 
result, and we use the termination exposing a (RU2O2) 
topmost layer as schematized in Fig. [21 as basis for our 
investigations. Previous work has shown that there are 
two prominent surface sites in each resulting rectangu- 
lar (lxl) surface unit-cell, which are primarily relevant 
for adsorption of O and CQ 6 i 7 i 8 i 9 i n i 13 i 14 i 15 i 16 i 17 i 18 i 19 i 20 . 
These are a bridge (br) site and a so-called coordinatively 
unsaturated (cus) site, the location of which is explained 
in Fig. [3 Depending on whether all bridge and/or all 
cus sites at the surface are either occupied by O or CO 
or empty, different (lxl) adsorption geometries result. 
For these we use a short-hand notation indicating first 
the occupancy of the bridge and then of the cus sites, 
e.g. O br /- for O adsorbed at all bridge sites, the cus 
sites being empty. 

Table [n] compiles the computed binding energies of O 
and CO at both surface sites in various of these con- 
figurations. The agreement between the full-potential 



LAPW/APW+lo and the standard PAW-std potential 
is at best modest. Although the potential seemed to 
perform satisfactorily in the standard tests, we now ob- 
serve deviations of up to 0.27 eV with respect to the full- 
potential binding energies. Interestingly, the PAW-/ po- 
tential that only marginally improved the / scattering 
properties performs clearly superior. Here, the agree- 
ment with the full-potential numbers is in all cases ex- 
cellent and the binding energetics is reproduced within 
±0.04 eV. On the other hand, the neglected relaxation 
of the semi-core states does not seem to play a big role. 
Relaxing the Ru Ap semi-core states as done in the PAW- 
Ap potential affects the binding energy values only little 
compared to the PAW-/ potential, cf. Table |nj A simi- 
lar result is obtained when further improving the O and 
C description through very hard small-core O and C po- 
tentials, as employed in the PAW-4p+ calculations. Also 
here, the binding energetics is little changed compared to 
the PAW-/ results. 

Investigating the reason for the severe error exhibited 
by the calculations employing the PAW-std potential, we 
first analyzed the density of valence states produced by 
the various PAW potentials and the LAPW/APW+lo 
method. Obtaining virtually indistinguishable results in 
all cases, we turned to the work functions as a measure of 
the average electrostatic potential. Table ITTT1 summarizes 
the obtained results, which neither provide a clear hint 
for the differing performance of the PAW-std potential 
compared to all others tested. While overall all PAW 
potentials represent the full-potential numbers quite well, 
there is in each case at least one configuration where 
we observe a deviation of the order of 0.1 eV. Changing 
from one PAW to the other, the description is sometimes 
improved and sometimes worsened, but this is not easily 
correlated with the binding energy changes listed in Table 

EH 

A correlation is, however, clearly identified when ana- 
lyzing the bonding geometries in more detail. The largest 
deviations between the full-potential and the PAW-std 
energetics occurs for the binding sites with the shortest 
Ru-0 or Ru-CO bondlengths. This concerns most no- 
tably the on-top adsorption of O cus at the cus sites, cf. 
Fig. [21 with a very short Ru-0 bondlength of 1.70 A in 
the O br /O cus phase. Also the Ru-O br bondlength is with 
1.92 A significantly shortened compared to the 1.99 A in 
RUO2 bulk, and the also badly described Ru-CO cus bond 
in the -/CO cus phase falls with 1.88 A in a similar range. 
These short distances then also reveal the origin for the 
different performance of the PAW-std and the PAW-/ 
potential: At such short distances, the tails of oxygen 
and carbon s and p states are picked up as high angular 
momentum contributions in the Ru PAW sphere. Since 
the Ru PAW-std potential is too repulsive for these chan- 
nels, cf. Fig. ^ the binding energy is underestimated at 
these sites. With this deficiency remedied in the PAW-/ 
potential through the differently chosen local component, 
the accuracy of the obtained surface energetics increases 
notably. This understanding also explains the seemingly 
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TABLE II: Binding energies (in eV/species) for and CO in (1 x 1) phases on RuO2(110) calculated with full-potential 
LAPW/APW+lo, and the various frozen-core PAW potentials (see text for the labeling). The binding energy of O atoms is 
calculated with respect to half the energy of a spin-polarized O2. 



Species 


Phase 


LAPW/APW+lo 


PAW-std 


PAW-/ 


PAW-4p 


PAW-4p+ 


Qbr 


O br /- 


-2.33 


-2.19 


-2.35 


-2.39 


-2.39 


Qbr 


O br /CO cus 


-2.16 


-2.05 


-2.18 


-2.20 


-2.20 


QCUS 


O br /O cus 


-0.86 


-0.59 


-0.85 


-0.87 


-0.89 


co cus 


o b vco cus 


-1.31 


-1.25 


-1.32 


-1.33 


-1.31 


co br 


co b v- 


-1.69 


-1.62 


-1.73 


-1.75 


-1.72 


co cus 


-/CO cus 


-1.48 


-1.38 


-1.50 


-1.52 


-1.50 



TABLE III: Work functions (in eV) of the (1 x 1) phases on RuO 2 (110) calculated with full-potential LAPW/APW+lo, and 
the various frozen-core PAW potentials (see text for labeling). 



Phase LAPW/APW+lo PAW-std PAW-/ PAW-4p PAW-4p+ 

-/- 3.83 3.86 3.85 3.83 3.83 

O b 7- 5.40 5.53 5.43 5.42 5.49 

O b 70 cus 7.01 7.05 7.05 7.03 7.00 

O b 7CO cus 5.83 5.96 5.96 5.97 5.92 

CO b 7- 5.02 5.06 5.06 5.06 5.02 

-/CO cus 5.29 5.30 5.31 5.31 5.27 



small differences between the PAW-std and the PAW-/ 
potential in the bulk tests. In the highly coordinated bulk 
Ru02 surrounding such short bondlengths never occur: 
They are specific to reactive, undercoordinated sites as 
typical for surfaces. The tiny imprecisions in the pseu- 
dopotential as reflected by the deviations in the scatter- 
ing properties of higher angular momenta did therefore 
not much affect the calculation of bulk properties, but 
have a large influence on the surface energetics. 



C. CO oxidation barriers at RuO2(110) 

These observations are, of course, not restricted to the 
binding of the probe molecules, but carry over to other 
surface quantities as well. To illustrate this point, we 
subject the reaction barriers for CO oxidation at the 
RuO2(110) surface to the same comparison with full- 
potential data. Adsorbed O and CO can react with each 
other if they occupy neighboring sites, which leads to four 
different possible reaction processes, namely O br +CO br , 
O br +CO cus , O cus +CO br and o cus +CO cus ^S. Fig. 03 
shows the initial state (IS) and transition state (TS) for 
two of these reactions. Since the reaction barrier results 
as the difference between the IS and TS energies, one 
would generally expect some degree of error cancella- 
tion, e.g. because the barrier is then no longer affected 
by any inadequacies in the atomization energies of the 
reacting species. From the data compiled in Table IIVI 
it is, however, obvious that the error due to the defi- 
cient local potential translates largely also to these quan- 
tities frequently targeted in catalytic studies. Again, the 
PAW-std potential underestimates all barriers by up to 
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FIG. 3: Perspective views of the RuO2(110) surface, show- 
ing the initial and the transition state configurations of the 
O b 7CO br and the O cu 7CO br reactions (Ru = light, large 
spheres, O = dark, medium spheres, C = light, small spheres). 



0.17eV, in particular those involving the shortly bonded 
qcus S p ec i es Q n the other hand, the agreement obtained 
between the PAW-/ and LAPW/APW+lo calculations is 
impressive, and is again not much affected through the 
test variations concerning relaxation of semi-core states 
(PAW-4p) or improved O and C description with harder 
potentials (PAW-4p+). Most crucial is therefore again 
the proper choice of the local potential, which if too 



7 



TABLE IV: Reaction energy barriers (in eV) for the four CO + O -^C02 reactions on the RuO2(110) surface, calculated in 
(1 x 2) unit-cells. Compared are full-potential LAPW/APW+lo calculations with the various frozen-core PAW calculations 
(see text for the labeling). 



Phase 


LAPW/APW+lo 


PAW-std 


PAW-/ 


PAW-4p 


PAW-4p+ 


o cus /co cus 


0.78 


0.66 


0.77 


0.78 


0.80 


o br /co br 


1.48 


1.46 


1.46 


1.47 


1.47 


o br /co cus 


0.99 


0.91 


0.97 


0.98 


0.99 


o cus /co br 


0.61 


0.44 


0.60 


0.61 


0.65 



repulsive as in the PAW-std case, leads to a too weak 
bonding of reaction intermediates with short bonds to 
the substrate atoms. Since these bonds are typically al- 
ready noticeably elongated at the TS configuration, this 
error is significantly reduced there. This way, cancella- 
tion can hardly occur in the difference between IS and 
TS, leading, for the PAW-std case, to reaction barriers 
that are too low. 

With this understanding we performed a series of cal- 
culations to assess the generality of our findings. There is 
no reason that our results should depend on the choice of 
the XC functional, and we indeed obtained the same un- 
dcrbinding of intermediates and underestimation of bar- 
riers when using another gradient-corrected (PW92) 35 
functional. On the other hand, a partial compensation of 
the error introduced by a PAW potential with a too re- 
pulsive local component could occur, when allowing the 
atomic positions to relax away from the equilibrium ge- 
ometry obtained within the LAPW/APW+lo approach. 
We tried this particularly for the most severely under- 
bound O cus species, but found only a negligibly increased 
binding energy even for the PAW-std potential. A sim- 
ilarly small effect was observed when all numbers were 
recomputed with another PAW potential where also the 
4s states were treated as valence. The variations within 
±0.02 eV with respect to the PAW-std calculations con- 
firm that in this system apparently neither the freezing of 
the 4s, nor of the 4p semi-core states has in the PAW ap- 
proach much influence on the accuracy of the computed 
surface energetics. The identified, much more critical in- 
fluence of the local potential on the other hand is unlikely 
specific to the PAW method. It will hold for other frozen- 
core approaches as well, and in particular for the closely- 
related ultrasoft pseudopotentials. In fact, we obtained 
results that were almost identical to the PAW numbers, 
when using USPPs constructed in a similar manner—, 
i.e. same local potentials and core radii. 

On this basis it would finally be nice to further ana- 
lyze the large discrepancies in hitherto published DFT 
numbers for this S y S tem£ii2ii£ii§iiLi2ii2i2!I. When doing 
so one has to recognize, however, that exploratory sur- 
face studies are normally not subject to such stringent 
convergence criteria as applied in this work. Insufficient 
basis sets could therefore equally account for the reported 
large scatter in the surface energetics. This is particu- 
larly a concern for the LAPW/APW+lo approach, where 
cutoffs as employed here rapidly lead to massive, if not 



prohibitive computational costs. Using cutoffs aiming 
only to obtain the correct energetic order in the reaction 
barriers, previous LAPW studies reported e.g. values 
that were up to 0.2 eV higher than the converged values 
found in the present worki&Si. In this respect we note 
in passing that in this system the two codes employed 
in this study converge from opposite ends with respect 
to the plane wave cutoff: in WIEN2k binding energies 
and reaction barriers become smaller with increasing ba- 
sis set size, while in VASP they become bigger. This 
leads rapidly to substantial differences in the absolute 
numbers at cutoffs aiming only to properly converge rel- 
ative energy differences within one approach. With this 
uncertainty with respect to the basis set convergence, 
we can at present only point out that our study shows 
that tiny deficiencies in the local component of the em- 
ployed pseudopotentials can also lead to substantial error 
in the surface energetics. This extra component could 
then possibly explain the particularly large differences 
in some published numbersi^i^i^ compared to the con- 
verged and perfectly agreeing full-potential and PAW or 
USPP values reported here. 



IV. SUMMARY AND CONCLUSION 

In conclusion, we have presented a detailed analy- 
sis of the accuracy of the frozen-core and pseudopo- 
tential approximation in DFT calculations at surfaces, 
by systematically comparing results obtained within the 
PAW method to those resulting from the full-potential 
LAPW/APW+lo approach. For the model system 
RuO2(H0) we find that the binding energies and reac- 
tion barriers of adsorbed O and CO are surprisingly lit- 
tle affected by the neglected relaxation of Ru semi-core 
states in the frozen-core PAW potentials. Instead, we 
identify a sensitive dependence of the surface energetic 
quantities on how well the choice of the local potential 
describes higher angular momentum scattering proper- 
ties. This is explained by the often very short bonds at 
undercoordinated surface sites, where the tails of adsor- 
bate s and p states reach noticeably into the substrate 
atom PAW sphere and are there picked up as high angular 
momentum contributions. If the local potential is then 
too repulsive, the binding is underestimated; if it is too 
attractive, the binding is overestimated. Since adsorbate- 
substrate bonds are often already substantially elongated 
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at transition state geometries, there is also little chance 
for cancellation of this error in calculated barriers. 

Highly-coordinated bulk environments exhibit rarely 
bonds that are as short as at reactive surface sites, so that 
corresponding small deficiencies in the local pseudopo- 
tential component are hardly noticeable, and might pass 
unnoticed in standard test calculations. In the present 
study, two choices of local potential led only to differences 
in the description of bulk energies by w 0.02 eV per O-Ru 
bond and to only seemingly irrelevant small differences 
in the scattering properties. On the other hand, they led 
to binding energies of shortly bonded O surface species 
that varied by 0.26 eV. We argue that this hitherto not 
much considered point may contribute largely to the of- 



ten substantially differing DFT results for this and other 
surface systems. 
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